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We present lattice calculations of kaon matrix elements with domain 
wall fermions. Using lattices with 6 /( 7 ^ = 5.85, 6.0, and 6.3, we estimate 
Bk{p ~ 2 GeV) = 0.628(47) in quenched QCD which is consistent with previ¬ 
ous calculations. At 6 / 5 ^ = 6.0 and 5.85 we find the ratio fx/'^p in agreement 
with the experimental value, within errors. These results support expecta¬ 
tions that 0{a) errors are exponentially suppressed in low energy {E a~^) 
observables, and indicate that domain wall fermions have good scaling behav¬ 
ior at relatively strong couplings. We also demonstrate that the axial current 
numerically satisfies the lattice analog of the usual continuum axial Ward 
identity. 


While lattice gauge theory has made signihcant progress in addressing the outstanding 
challenge of calculating hadronic observables from first principles, a basic feature of the 
strong interactions has been missing in these calculations, the SU(Nf)L x SU(Nf)R chiral 
flavor symmetry of the light quarks which is broken explicitly by present lattice discretiza¬ 
tions of continuum QCD. We recently reported [jl[ on calculations using a new discretization 
for simulations of QCD, domain wall fermions (DWF) [^,Q, which preserve chiral symmetry 
on the lattice in the limit of an infinite extra 5th dimension. There it was demonstrated 
that DWF exhibit remarkable chiral behavior |l[] even at relatively large lattice spacing and 
modest extent of the fifth dimension. Here we give further results using DWF which are of 
direct phenomenological interest. 

In addition to retaining chiral symmetry, DWF are also “improved” in another important 
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way. In the limit that the number of sites in the extra dimension, Ng, goes to inhnity, the 
leading discretization error in the effective four dimensional action for the light degrees of 
freedom goes like O(a^), unlike the case for ordinary Wilson fermions, for which the errors 
are 0(a), a being the lattice spacing. This theoretical dependence is deduced from the fact 
that the only operators available to cancel 0(a) errors in the effective action are not chirally 
symmetric; thus no 0(a) errors exist in the low energy theory. For hnite Ng, 0(a) corrections 
are expected to be exponentially suppressed with the size of the extra hfth dimension. Our 
calculations for show a weak dependence on a that is easily £t to an ansatz. At 
(3 = 6.0 (/3 = Q /the lattice spacings determined from rrip and /^r agree within less than 
five percent. This improved scaling behavior is plausible in light of the fact that DWF retain 
an important continuum symmetry at non-zero lattice spacing. 

As in our previous paper, we use the boundary fermion variant of DWF 0] developed by 
Shamir |^. The DWF action is essentially a hve dimensional analog of the ordinary Wilson 
fermion action with two key differences: (1) the relative sign between the Wilson term and 
the (hve dimensional) Dirac mass, M, is opposite to the usual convention. This leads to the 
appearance of massless chiral modes on the boundaries of the fifth dimension, a left handed 
fermion on one wall and a right handed one on the other. (2) The layers s = 0 and s = Ng — 1 
(s denotes the coordinate in the extra dimension) are coupled with strength —m (m > 0). 
Neglecting exponentially small corrections, in Ref. 0 it was shown that the parameter m 
is (proportional to) the mass of the light four dimensional quark which is assembled from 
the two chiral modes, niq = mM(2 — M). The chiral limit is A^* oo and m — >■ 0, which 
requires no hne tuning unlike ordinary Wilson fermions. 

Recently the exponentially small corrections to the quark mass have been given at tree 
level 0, TUq = M{2 — M)(m + {I — M + 0(p‘^)}^‘’). In the presence of interactions M 
is renormalized additively just like ordinary Wilson fermions which also acquire a “mass 
term” proportional to Perturbatively at one loop the main effect of the interactions is 
the replacement M ^ M + g‘^'E(p) where S(p) is the quark self-energy. Thus the chiral 
limit still holds with the replacement (1 — > (Me — . This is analogous to the 
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renormalization of the critical hopping parameter from its tree level value of 1/8 in the case 
of ordinary Wilson fermions. Of course, the crucial difference is that for DWF the additive 
corrections are exponentially suppressed. In our original study we found Me ~ 1.7 for non- 
perturbative couplings corresponding to quenched simulations at /? ~ 6.0, which also agrees 
roughly with a simple mean field argument [Q. 

For QCD, the DWF are gauged in the ordinary four dimensions only, and the left and 
right handed modes couple equally to the gauge field. Thus the five dimensional theory 
gives rise to a low energy effective theory {E « a“^) describing interacting vector quarks 
in four dimensions whose right and left handed components are localized around s = 0 and 
s = Ng — respectively. 

In Ref. P] it was shown that operators constructed from the quark helds formed by the 
chiral modes on each wall satisfy the following four dimensional chiral Ward identities(CWI). 

2/2, •••)) = 2m (J5 (a;)0(2/i, 2/2,...)) 

+ ^(<^10(2/1,2/2, ■••)), ( 1 ) 

which result from demanding invariance of (0(2/1, 2/2, •••)) under an inhnitesimal axial trans¬ 
formation Here is the four dimensional partially conserved axial current (PCAC) 
constructed from a sum of hve dimensional helds over all slices in the extra dimension [^. 
The operators 0(2/1, 2/2, •••) and the pseudoscalar density J^{x) are constructed from four 
dimensional quark helds using the chiral modes on each boundary, 

^ 0) -1- ^ Ns - 1), 

q{x) ='ip{x,Ns-iy^-^+ (2) 

In Eq. ^ Jgg is an anomalous pseudoscalar density that results from the non-invariance of 
the action under the axial transformation for hnite Ng. For havor non-singlet currents, this 
contribution to the r.h.s. of Eq. vanishes identically in the limit W —^ C)0 |^, and we 
are left with the continuum-like relations. Below we demonstrate explicitly that aX j3 = 6.0 
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Eq. 2 is satisfied for the usual PCAC Ward identity, 0 {yi,y 2 , ■■■) = J 5 , and the anomalous 
contribution is small for Ng = 10 and reduces further by more than a factor of 2.5 as Ng is 
increased to 14. 

To obtain Bk we need the matrix element of the As = 2 four quark operator that governs 
K — K mixing, 0{y)LL = [s7j/(1 — 75 )^]^, which depends on the energy scale y. On the 
lattice and using DWF, Oll is a simple transcription of the above using the quark fields 
in Eq. |^. Sandwiching Oll between K and K states and taking the ratio with its value in 
vacuum saturation yields Bk- 

TABLE I. Summary of simulation parameters, “size” is the number of spatial sites times the 
temporal extent times Ng. M is the five dimensional Dirac fermion mass, and m is the coupling 
between layers s = 0 and — 1. The number in parentheses is the number of configurations used 
at each value of m. 


(3 

size 

M 

m(^ conf) 

5.85 

16^ X 32 X 14 

1.7 

0.075(34), 0.05(24) 

6.0 

16^ X 32 X 10 

1.7 

0.075(36), 0.05(39), 0.025 (17) 

6.3 

243 X 60 X 10 

1.5 

0.075(11), 0.05(14) 


In order to investigate the continuum limit of quenched QCD with DWF, we have carried 
out simulations at gauge couplings [3 = 5.85, 6.0, and 6.3. The simulation parameters are 
summarized in Table |. The number of configurations in our study is rather small, and we 
have made no attempt to estimate finite (four dimensional) volume systematic errors. These 
deficiencies will, of course, be addressed in future works. The lattices correspond to (1.5 fm)^ 
for (3 = 6.0 and 6.3 and (2.1 fm)^ for 5.85, which from previous lattice studies may result 
in a deviation of a few percent from the infinite volume case. We have begun to address 
systematic corrections due to finite Ng. All of the correlation functions discussed below were 
calculated in the lattice Landau gauge which was chosen for convenience, in principle any 
gauge will do. 

We begin by discussing the numerical investigation of Eq. for O = J5. First, Eq. |l| 
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is satisfied exactly on any configuration since it is derived from the corresponding operator 
identity. We checked this explicity in our simulations. In the asymptotic large time limit we 
get 


2 sinh(am^/2)(A^|7r)/(J5|7r) = 2m + 2(J5g|7r)/(J5|7r), (3) 

which goes over to the continuum relation for am^r 1 and Ng ^ oo. At /3 = 6.0 and 
Ng = 10 we hnd the l.h.s. of Eq. || to be 0.1578(2) and 0.1083(3) for m = 0.075 and 
0.05, respectively. The anomalous contributions for these two masses are 2x(0.00385(5) 
and 0.00408(12)), which appears to be roughly constant with m. Increasing Ng to 14 at 
m = 0.05, the anomalous contribution falls to (2x) 0.00152(8) while the l.h.s. is 0.1026(6), 
which shows that increasing Ng really does take us towards the chiral limit. 

While we have not investigated the CWI for Oll^ fhe matrix element (^K^\Oll\K^^ 
vanishes linearly with m in the chiral limit as required by chiral perturbation theory and 
shown in Fig. |^. This indicates that the anomalous term in Eq. |^for Oll is highly suppressed. 
At j3 = 5.85 the two data points extrapolate linearly to -0.0005(100) at m = 0. At (3 = 6.0 
the three data points extrapolate to -.004(9) with a degree of freedom of 0.2 for a 

non-covariant £t. We do not have enough data to perform a covariant fit. At (3 = 6.3, the 
two points extrapolate to 0.05(3). This slight overshoot is not unexpected since the values 
of m used are for rather heavy quarks. In our initial study we found a similar behavior [|I|, 
and as the quark mass was lowered, the required linear behavior set in. Since we do not 
have a smaller mass at this coupling, there will be a small systematic increase in since 
the fit will overestimate the matrix element at the strange quark mass. All of the above 
results are for Ng = 10 except ai j3 = 5.85 where Ng = 14 was used for reasons explained 
below. 

Fig. 1^ shows Bk as a function of amp, or equivalently the lattice spacing. B^ is estimated 
at each coupling from a linear £t of the degenerate quark data. The £t is then evaluated at 
one-half the value of m corresponding to the strange quark as determined from a fit to the 
pseudoscalar mass squared (see below), arup is determined from a simple jackknife average 
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of the effective mass over a suitable plateau. The results for Bk depend weakly on (3, and 
are well £t to a pure quadratic in a. Using this fit, we hud = 0.628(47) in 

the continuum limit with a conhdence level of 0.39. This result is already consistent with 
the previous Kogut-Susskind result [^,0 and a recent Wilson quark result ^ using CWI’s 
similar to Eq. to enforce the proper chiral behavior of Oll for Wilson quarks. We note 
that the data can be £t to a linear function of the lattice spacing as well, which yields 
Bk{h = = 0.617(80), though we emphasize again that linear corrections are expected 

to be highly suppressed on theoretical grounds. More importantly, there is no evidence for 
0(a) corrections in Fig. |I|. Similarly, the denominator in Bk-, (K\Oll\K\ , exhibits the 
correct chiral behavor. Above, the notation = a“^) simply means the uncorrected 

lattice data have been used to perform the extrapolation; i.e. , our result does not include 
the perturbative running of B^ at each lattice spacing to a common energy scale. This 
requires a perturbative calculation to determine the renormalization of Oll, which has not 
yet been done. The energy scale at /? = 6.0, as determined by the inverse lattice spacing, is 
roughly 2 GeV. Since the dependence on a is already mild, this should not have a signihcant 
impact on our a = 0 estimate. We also note that the coefficient of the quadratic term in our 
fit is 0.12(23), or zero within errors. Also, from our previous work |]^ which did not give a 
value for Bk, we find BxifJ^ = a“^) = 0.67(4) on a set of 20 Kogut-Susskind lattices with 
"^Ks = 0.01 and (3 = 5.7 0 and the same five dimensional lattice volume as the point at 
(3 = 6.0. The energy scale is nearly that oi (3 = 6.0, quenched. This partially unquenched 
result indicates that the error from quenching may be small, as was found for Bk using 
Kogut-Susskind quarks |T0|JTT[ . 


At [3 = 6.0, we have also calculated Bk using the partially conserved axial current A^{x) 
(and the analogous vector current) at m = 0.05 and 0.075. This point split conserved 
current requires explicit factors of the gauge links to be gauge invariant. Alternatively 
a gauge non-invariant operator may be defined by omitting the links; the two definitions 
become equivalent in the continuum limit. Results for the gauge non-invariant operators 
agree within small statistical errors with those for operators constructed from Eq. H. The 
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results for the gauge invariant operators are somewhat larger: = a~^) = 0.872(22) 

and 0.926(19) at m = 0.05 and 0.075, respectively. A similar situation holds in the Kogut- 
Susskind case where it was shown that the gauge invariant operators receive appreciable 


perturbative corrections which bring the two results into agreement |10 


Using Eq. |^, neglecting the anomalous contribution, and using the dehnition of the 
decay constant fpmp = (0|Aq|P), we can determine the pseudoscalar decay constant from 
the measurement of (0| Jg |P). Performing simultaneous covariant hts to the wall-point and 
wall-wall correlators of J 5 yields the matrix element. At /3 = 5.85 and 6.0 each £t has a 
good conhdence level (CL 0.7). The results are shown in Fig. Proceeding as before 
with Bk, we hnd fx = 159(14) MeV and 164(12) MeV for (3 = 6.0 and 5.85, respectively. 
The errors are statistical and do not include the error in the lattice spacing determination 
from arUp. The central values agree with experiment, fK+ = 160 MeV. The lattice spacing 
determinations from arUp give a~^ = 1.53(27), 2.09(21), and 3.20(81) GeV at /3 = 5.85, 6.0, 
and 6.3, respectively. These are similar to Wilson and Kogut-Susskind lattice spacings for 
similar quenched lattices. Alternatively, we may form the dimensionless ratio We 

hnd for (3 = 5.85 and 6.0, fx/^p = 0.213(42) and 0.206(27), where we have added the 
statistical errors naively in quadrature. The experimental result is 0.208. At present our 
data at /3 = 6.3 are too noisy after extrapolation to give a significant result. Finally we 
note that we have also calculated the decay constant directly from the matrix element of 
the partially conserved axial current at the points m = 0.05 and 0.075 aX j3 = 6.0, and the 
results agree with those using the matrix element of the pseudoscalar density (see Fig. 
While the above results indicate good scaling behavior, they must be checked further with 
improved statistics and a fully covariant htting procedure. More importantly, the continuum 
limit still has to be taken: a recent precise calculation using quenched Wilson quarks by the 
CP-PACS collaboration gives a value for fx/^p in the continuum limit that is inconsistent 


with experiment [|^], so the above agreement with experiment may be fortuitous. 

In Fig. ^ we show the pion mass squared as a function of m. Lowest order chiral per¬ 
turbation theory requires to vanish linearly with m. At /3 = 6.0 our three data points 
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extrapolate linearly to .008(10), which agrees with the above expectations. At /3 = 6.3, the 
two data points extrapolate to —0.021(3), which is once again likely due to quark masses 
that are too heavy to agree with lowest order chiral perturbation theory. At /d = 5.85 the 
two masses extrapolate to 0.045(10) for Ng = 10 and 0.031(13) for Ng = 14. This discrep¬ 
ancy is probably not due to higher order terms in the chiral expansion since the physical 
quark masses are light compared to the masses at the other couplings, and the curvature 
would have the wrong sign. We see a large downward shift in as Ng goes from 10 to 
14. However, increasing Ng to 18 at m = 0.075 has a negligible effect. This behavior may 
signal a strong coupling effect where the suppression of explicit chiral symmetry breaking 
terms with Ng may be weakened. In the case of the vector Schwinger model, it was found 
that topology changing gauge conhgurations can induce signihcant explicit chiral symmetry 
breaking effects [^]. Further investigation is required. 

Our study shows that DWF are an attractive alternative to Kogut-Susskind and Wil¬ 
son quarks for lattice QCD calculations where chiral symmetry is crucial. For weak matrix 
elements in particular, DWF yield good agreement with expectations from chiral pertur¬ 
bation theory without the complicated mixing of operators required with Wilson quarks, 
or the entanglement of flavor and space-time degrees of freedom as with Kogut-Susskind 
quarks. Perhaps even more importantly, up to exponentially small corrections, DWF main¬ 
tain the full chiral symmetry of QCD at relatively strong couplings, and thus should have 
more continuum-like behavior. The data presented here seem to indicate just that, though 
future studies with improved statistics are needed to conhrm this. This improved scaling 
may compensate for the added cost of the extra dimension. 


Our domain wall fermion code relies heavily on the four dimensional MILC code [13 


which we are happy to acknowledge again. This research was supported by US DOE grant 
DE-AC0276CH0016. The numerical computations were carried out on the NERSC T3E. 
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FIG. 1. The matrix element {P\Oll\P) vs. m. \P) is a non-singlet pseudoscalar state, m 
is proportional to the quark mass in lattice units. At /3 = 5.85 and 6.0 (P\Oll\P) vanishes 
linearly with m. The slight overshoot at /? = 6.3 is likely due to higher order terms in the chiral 
expansion(see text). 
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FIG. 2. Kaon B parameter. The solid line is a pure quadratic ht to the data, and the burst 
denotes the extrapolation to the continuum limit, a = 0. The data are for Ng = 10 (/? = 6.0, 6.3) 
and 14 (5.85). The cross (not used in the fit) denotes the partially unquenched result discussed 
in the text. The energy scale is roughly 2 GeV at /3 = 6.0, and perturbative corrections have not 


been included. 





FIG. 3. Pseudoscalar decay constant. Bursts are linear extrapolations to values of m corre¬ 
sponding to the pion and the kaon. Ng = 10 and 14 for /3 = 6.0 and 5.85, respectively. The crosses 
denote values calculated from the matrix element of the partially conserved axial current (/3 = 6.0 




FIG. 4. Pion mass squared. The curve at (5 = 6.0 {Ng = lO(octagons)) agrees with the expecta¬ 
tion from lowest order chiral perturbation theory. The others show small but significant deviations. 
At P = 5.85 going from Ng = 10 to 14 reduces the discrepancy. However increasing Ng to 18(cross) 
has no effect. At P = 6.3, the discrepancy is most likely due to higher order terms in the chiral 


expansion (see text). 




